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^D Abstract. We determine numerically the self-similarity maps for spectra of the 

\^ almost Mathicu operators, a two-dimensional fractal-like structure known as the 

— 1 Hofstadter butterfly. The similarity maps each have a horizontal component deter- 

mined by certain algebraic maps, and vertical component determined by a Mobius 
transformation, indexed by a semigroup of the matrix group G'L2(Z). Based on the 
numerical evidence, we state and prove a continuity result for the similarity maps. 
^H We note a connection between the indexing of the similarity maps and Morita 

f — equivalence of rotation algebras Ag, a continuous field of C*algebras. 
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Figure 1. The Hofstadter butterfly, a fractal-like structure of line spectra. 



1. Introduction 

Looking at the image presented in Figure 1, one immediately notices the striking, 
repetitive pattern of "butterfly wings" that march off towards the vertical horizons at 
the top and bottom of the image. This rendering of the so-called Hofstadter butterfly, 
drawn to high-resolution using a combination of numerical algorithms and PostScript 
graphics programming, reveals some beautiful symmetries of a fundamentally math- 
ematical object. The goal of this paper is to specify exactly the symmetries of this 
image and prove continuity results of the corresponding similarity maps, motivated 
by numerical evidence collected in our study of the butterfly. 
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Although the Hofstadter butterfly looks like a fractal, the image is not created 
using the usual iterative equations or recursive graphical methods designed to visually 
render typical fractals (Mandelbrot, 1990). Rather, this image is created from an 



explicit numerical computation of spectral values of a family of linear operators on 
Hilbert space. 

More precisely, this image is constructed as a layered suite of horizontal lines 
which arise as the spectra of almost Mathieu operators. Each of these operators is 
succinctly represented as the self-adjoint element 



(1) 



hf, 



u + u + V + V 



in the C*-algebra Ag generated by universal unitaries u, v satisfying the commutation 
rule 



(2) 



vu 



e^^'^uv. 



In Figure 1, the vertical axis is spanned by the parameter 9, with Q < 9 < 1, while 
the horizontal axis corresponds to spectral values x in the range —4 < x < 4. 

Only rational values of 9 are used in the construction of the image in Figure 1. For 
this reason, the image is properly called the rational Hofstadter butterfly. Extensions 
to irrational values of 9 has been a theme in the long history of study these operators 
and the physical problems that motivated it, going back at least to a mathematical 
analysis of Bloch electrons (Brown, 1964). A selection of relevant studies over the 



years is given in the references, including (Avila and Jitomirskaya, 2006; Bellissard 



and Simon, 1982 Bellissard, 1990; Choi et al., 1990 Connes, 1994 Goldman, 2009 



Hofstadter, 1976 Kaminker and Putnam, 2003; Last, 1994 Puig, 2004 Ypma, 2007) 



Unlike the irrational case as considered in (Arveson, 1994), rational values of 



9 = p/q lead to the computation of spectral lines based on finite dimensional eigen- 
value calculations using tridiagonal qx q matrices (Lamoureux, 1997). In particular. 



the endpoints of the spectral lines are specified by eigenvalues of these matrices. 
Numerical algorithms for the tridiagonal eigenvalue computations are rapid and ac- 
curate, making calculations of the rational butterfly ideal for numerical experiments. 
Based on these experiments, we deduce a systematic catalogue of the symmetries of 
the butterfly. 

Repetition of a geometric object fading out to infinity is characteristic of hyperbolic 
geometry (Coxeter, 1942), hence it is perhaps no surprise that the matrix group 
GL2(Z) arises in the observed symmetries. Sections 2 through 5 reveal a semigroup 
of linear fractional transformations represented by elements of GL2(Z) which act on 
the vertical parameter 9 to generate the family of similarity maps. 

Equally important are the algebraic curves specifying the horizontal component of 
the similarities, acting on the horizontal parameter x. Sections 6 through 9 present 
numerical evidence of the continuous maps taking one level of horizontal spectra 
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to another, using an indexing of intervals and a correspondence of image points 
under polynomial maps. These are precisely the characteristic polynomials of the 
aforementioned finite dimensional g x g matrices giving the spectra. 

Uniting the vertical and horizontal components leads to an identification of the 
general form of the similarity maps, as presented in Section 10. The similarity extends 
to gap-labelling, a specific method for indexing the butterfly wings based on Chern 
characters, which is discussed in Section 11. Proof of continuity of the similarity maps 
is given in Section 12 and generators for the semigroup of similarities is presented in 
Section 13. 

"Xenocides, who is ugly, makes ugly poetry," said Aristophanes. With the Hofs- 
tadter butterfly, we see pretty mathematics making pretty pictures. 

2. Numerical evidence: Similarity maps 

Mirror image symmetry in the horizontal direction suggests an obvious self-similarity 
of the image in Figure 1 given by a reflection about the line x = 0, namely: 

(3) {x,e)^{-x,e). 

This map is continuous, and does map the butterfly to itself, as the spectrum of each 
hg is symmetric. 

In the vertical direction, another mirror image symmetry in the figure suggests a 
second self-similarity map given by reflection about the line 6 = 1/2, namely: 

(4) {x,e)^{x,l-0)- 

Again, this is a continuous map, and again maps spectra properly since the algebra Ag 
is isomorphic to Ai_g, with operator hg mapping onto hi_g under the isomorphism. 
Now, a more interesting symmetry is observed mapping the large central butterfly 
onto the next largest butterfly in the bottom half of the image in Figure 1. Zooming 
in on this butterfly, as shown in Figure 2, one sees the top of the butterfly at ^ = 1/3, 
the bottom at 6' = 0, and the centre of the butterfly at 6^ = 1/4. This suggest we 
must find a self-similarity map that, on vertical parameter values 6, will map 9 ^^ 9' 
as 

(5) h-> 0, 1/2 ^ 1/4, 1 ^ 1/3. 

No linear map will do, but a linear fractional transformatioEJj does work, using the 
map 

(6) 9^9' ^ 



29 + 1 



i.e. a Mobius transformation 
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Figure 2. Similarity to the bottom third of the butterfly. 



It is reassuring to notice that other obvious lines in the full image map correctly 
under this map. For instance, the line at ^ = 1/3 in the full butterfly in Figure 1 
should map to the top of the second largest butterfly in Figure 2, which is at ^' = 1/5. 
And indeed, the LFT here does just that. Similarly we can check for that the line 
at 1/4 maps to 1/6, the line at 1/5 maps to 1/7, and so on. Checking numerically 
many of these lines assures us that the linear fractional transformation is the proper 
choice. 

Going to a numerical experiment based on these observations, we construct a map 
from the full image in Figure 1 to subset like Figure 2, using the LFT in the vertical 
component, and contracting linearly in the horizontal component. Specifically, we 
map according to the rule 



(7) 



(x, 6) ^ (x', e') 



;i-.82*0)x. 



9 



26 + 1 



The coefficient 0.82 was chosen to get the correct width for the top horizontal line. 
The result is the image shown in Figure 3, which is very much like the lower butterfly 
image in Figure 2. 

Observing the slight differences between the image in Figure 1 and the result of 
the numerical experiment shown in Figure 3, (the difference in the curvature of the 
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Figure 3. A rendering of an approximate similarity to the bottom 
third, using a hnear fractional transformation vertically, linear scaling 
horizontally. 

wingtips for instance), we can conclude the whatever the similarity map is, it is 
only approximately linear. Nevertheless, the basic structure and position of the tiny 
butterflies are preserved. This serves as numerical evidence that the vertical maps 
are indeed given by linear fractional transformations. The horizontal component is 
more complicated, and is discussed in Sections 7, 8, 9. 

3. Numerical evidence: Vertical similarities from GL2{Z) 

Perhaps it is suggestive that two of the self-similarity maps identified so far (verti- 
cal flip, and the bottom third map) have vertical components give by Mobius trans- 
formations, 

1 n n 

(8) e ^ and 9 ^ 



1 29 + 1 

Both transformations are specified by a 2 x 2 matrix in G'L2(Z), in the form 



(9) M 



a b 
c d 



yielding a corresponding linear fractional transformation 
(10) 6^9' "^ + ^ 



c9 + d' 

It will be convenient to label these maps by the corresponding matrix, and we note 
that composition of the LFTs corresponds to matrix multiplication in G'L2(Z). 

A leap of faith suggests looking for symmetries in the butterfly indexed by elements 
of GL2(Z). Examining Figure 4, we see the bottom of the central core of butterflies, 
and we can compute a list of matrices that implement the self-similarity map on the 
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Figure 4. A sequence of similarities in the lower central core of the 
butterfly, and matrices in SL^i^I^ that implement the vertical compo- 
nent of the similarity. 



vertical component. For instance, we can map the whole butterfly onto the central 
bottom butterfly that tops out at label Q = 1/5. For this, we need an LFT that 
maps H- 0, 1/2 i-> 1/6, 1 H- 1/5. Solving for coefficients a,b,c,d in the LFT, we 
find the matrix 



(11) 



1 

4 1 



will work to implement an LFT that maps the whole figure into this central bottom 
butterfiy. 

Continuing down the central core, we see an infinite sequence of butterfiies extend- 
ing to the "horizon" at ^ = 0. Some simple calculations analogous to those above 
produce a sequence of similarity maps indexed by matrices of the form 



(12) 



1 
2n 





1 





1 



2ra 



A"", for n > 1. 



We will see later that this matrix A is a key generator of the similarities. 

Hopping to the top of butterfiy, in Figure 5 we see a central core of butterfiies 
extending to the upper "horizon" at level 6 = 1. Here, the self-similarity maps are 
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Figure 5. A sequence of similarities in the upper central core of the 
butterfly, and matrices in GL2(Z) that implement the vertical compo- 
nent of the similarity. 



obtained using matrices in GL2{Z) of the form 



(13) 



l-2n In 
-2n 2n + l 



for n > 1. 



Matrices from symmetries on top (Eqn 13) are related to matrices for the symme- 
tries on the bottom (Eqn 12), via the conjugation 



(14) 



1 -2n 
-2n 



1 

2n + l 



= 


'-11" 

1 




' 1 " 

2n 1 




'-11" 
1 



This is just conjugation with the flip symmetry map of the butterfly, as the matrix 
(15) B = 



-1 1 
1 



induces the linear fractional transformation mapping 6 to 1 — 6, turning the butterfly 
upside down. 

Attending to some of the butterflies on the side of the image, we flnd the relevant 
linear fractional transformations are represented by matrices of the form 



(16) 



1 
-1 n 



for n < 3, 
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Figure 6. A sequence of similarities on the edge of the butterfly, and 
matrices in GL2{'L) that implement the similarity. 



which generate the symmetries on the bottom half of the image shown in Figure 6. 
On the top of the image, we obtain transformations indexed by matrices of the form 



(17) 



n 
n 



2 
1 



for n > 3. 



Again, we observe there is a relation between these two forms of transformation, 
through conjugation with the flip, as we have 



(18) 



1 n-2 
1 n- 1 



-1 1 
1 



1 
-1 n 



-1 1 
1 



10 
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Note, however, that these self-similarity maps have a certain discontinuity. In 
particular, the centre of the large butterfly gets mapped to a "broken" butterfly on 
these side maps. This tells us the self-symmetry maps need not be continuous. In 
particular there may be a discontinuity in the horizontal direction as we cross the 
X = spectral value. 

Nevertheless, this is a mild discontinuity and is easy to understand. In the con- 
struction of the Hofstadter butterfly, for 9 = p/q with q even, the spectrum consists 



of q intervals on the real line, two of which touch at the point zero. (Choi et al., 1990) 



It is these "touching intervals" that are getting broken apart in the above similarity 
maps. So although there is an apparent discontinuity, it is perhaps better to describe 
it as the breaking apart of two touching spectral lines. The similarity map should 
take the single, common endpoint, and map it to two distinct endpoints of two non- 
overlapping intervals. The details of this double-valued map will be discussed in 
Section 10. 



4. Numerical evidence: A symmetry with break 

As noted in the previous section, under the similarity map sometimes the butterfly 
breaks. This suggests we can look for more similarities if we are a bit more open to 
what a broken butterfly looks like. 

In Figure 7, we have an example of a really broken butterfly. The butterfly fits 
into the region 

(19) 0<9< 1/2, 

and the similarity map is given by the linear fractional transformation with matrix 



(20) 



1 

1 1 



This break (at 6 = 1/3) looks pretty bad, but in fact we know from gap labelling 
theorems that when one plots the gaps an inverse slope 2, there is a discontinuity 
exactly at 6' = 1/3 (Lamoureux, 2010). So in fact this is just confirming the fact 



that the butterflies are coming from gap labelling. 
By including this symmetry, we include the matrix 



(21) 



A 



1 
1 1 



in our collection of matrices in (7^2 (Z) generating similarity maps. 

We observe that all the LFTs seen so far map the 6 interval [0, 1] into itself, and 
are indexed by certain elements of of GL2(Z). This evidence suggests a certain 
semigroup in GL2(Z) specifies the possible self-similarity maps. A precise statement 
identifying the semigroup is given in the next section. 
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Figure 7. An unusual, cracked similarity. Note the basic butterfly 
structure spanning the bottom half of the full image. 



5. Generating the GL2{Z) symmetries 

The following technical result states that the set of linear fractional transforma- 
tions which map the interval [0, 1] into itself forms a semigroup, indexed by a 2- 
generated semigroup of matrices in GL2{'L). 

Theorem 1. Let G he the semigroup of linear fractional transformations of the form 

aO + h 



(22) 



e^e' 



c9 + d' 



for some 



a b 
c d 



G GL2(Z)/{±/}, 



which map the interval [0, 1] into itself. Then G is generated by the two maps 
6/{6 + 1) and ip{6) = 1 — 9, represented by matrices 



(23) 



A 



1 
1 1 



and B 



-1 




The proof of this technical result is given in Appendix 1, and is based on the 
Euclidean algoritm. We make a few observations about this theorem. 

The LFTs do not notice a change of sign in the representing matrix, so the theorem 
refers to the quotient group G'L2(Z)/{±/}. 



12 
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Figure 8. Multiple similarities, indexed by the same matrix in SL2{'L) 



These two generators correspond to two symmetries for the Hofstadter butterfly. 
The first generates the map that takes the butterfiy to its lower half, mentioned in 
Section 4. The second is the vertical fiip symmetry, mentioned in Section 2. By 
including the vertical fiip in our symmetries, we are able to generate all (vertical 
components) of the symmetries using only two generating maps. 

The fact that the Euclidean algorithm is involved suggests that determining the 
details of a symmetry map will be rather involved - the factorization into generators 
is not simple. 

It is worth mentioning a deep result from C*-algebras, that states two rotation 
algebras Aq, Agi are Morita equivalent if and only if there is a linear fractional trans- 
formation mapping 



(24) 



e^e' 



aO + h 
cO + d' 



for some matrix 



a 
c 



b 
d 



connection this might have with the 



in GL2{'L) (Effros and Shen, 1980). We have no idea what 



self-similarity maps we have above, which act 



on spectra of operators in the algebras Aq^ Aqi, not on the algebras themselves. 



6. Horizontal similarity: Interval maps 
Considering now the horizontal component of the similarity maps, a close exami- 



nation of the apparent symmetries suggests the intervals /i, I2 



, Iq in the spectrum 



at level 6 = p/q are mapped onto a subcollection of the intervals I[, I2, . . . I' , in the 
spectrum at level 6' = p' /q'. We note here which sub collections appear. 
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Examining Figure 8, notice for the single vertical map with linear fractional trans- 
formation given by matrix 



(25) 



there are in fact three choices of horizontal maps. Start first at initial 6 = 1, which 
maps to 6' = 2/5 under this LFT. The spectrum at level 6 = 1 has one interval, 
while at level 6' = 2/5, there are five intervals that we can map to. Only three are 
obtained. In the left part of the image in Figure 8, the initial interval Ji maps onto 
interval I[; the middle image maps onto interval I'^, and the right image maps onto 
interval Jg. Writing p' /q' = 2/5, we see that the three interval maps are given as 



T' 

-'Op'+l' 



h^l'= /' 



Ip'+l) 



h ^ I', 



T' 

^2p'- 



(26) h^I[ 

Checking now at the level 6 = 1/2 mapping to 6' = 3/8, the two intervals /i, I2 at 
6 = 1/2 map to two interval I[, I2 on the left; to the two interval I'^, I^ in the middle; 
and Ij, I'g on the right. In this case, we have p' /q' = 3/8 and we can summarize the 
left, middle and right maps as 

(27) Ik I— ^ I^pi+ki ^k ^^ hp>+ki ^k I— ^ hp'+ki 

for k = 1,2. Similarly, at level ^ = 0, the single interval Ji maps as 

(28) 1,^/; = /^^,^^, h^r, = i[^,_^„ 1,^/^ = 4,^^, 

where p' /q' = 1/3. The form of the interval map is very consistent. 

We go back and check the other similarity maps. In Figure 4, the interval maps 
are apparently of the form 

(29) Ik ^ C+k = 

for the vertical similarity given by 

1 

2n 1 



ji 

n-p' +ki 



(30) 



We can see this by observing the q intervals at level 6 = p/q map to the central q 
intervals of the set of q' = 2np + q intervals at level 6' = p' /q. 
In Figure 5, the interval maps are of the form 



(31) Ik *-> In-{q-p)+k 

for the vertical similarity given by 

1 -2n 



/: 



n-{q'—p')+k' 



(32) 



-2n 



2n 
2n + l 



Again, the q intervals at level 6 = p/q map to the central q intervals of the set of 
q' = 2n{q — p) + q intervals at level 6' = p' /q. 
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In Figure 6, a similarity in the top half, say with matrix 

(33) 



1 1 
1 2 



gives interval maps of the form 

(34) 4 ^ /;y+fe, 

with r = 1 on the right side of the butterfly (as illustrated), and r = on the left 
side of the butterfly (not shown in Figure 6). However, a similarity in the bottom 
half of Figure 6, such as with matrix 



(35) 



1 
-1 3 



gives interval maps of the form 

(36) Ik ^-^ Ir-(q'-p')+ky 

with r = 1 on the right map, r = on the left map. 

Generally speaking, it appears there are two types of interval maps, those of the 
form 



(37) 


4 ^ I'r.p>+k, 


and those of the form 




(38) 


J-k ^^ ir-{q'-p')+ki 



for some choices of non-negative integer r. 

In the next three sections we determine how the horizontal component of the 
similarity map behaves within each interval 4. 

7. Horizontal similarity: Cubic case 1 1-> 1/3 

The vertical component of the self-similarity map is given by linear fractional 
transformations, as described in Sections 2 through 4. The horizontal component 
maps intervals to intervals in the spectra, and appears to be nearly linear, as we saw 
in the construction of Figure 3. We investigate the details of this horizontal map. 

Again referring to the butterfly list in Figure 1, we see that the spectral line at 
height ^ = 1 is approached by a dense cluster of spectral points a heights 9 < 1. 
Under the LFT mapping 6 \-^ 6' = 6/ {26 + 1), the point 6 = 1 maps to 6' = 1/3, and 
this dense cluster of spectral points near 1 maps to a dense cluster near 1/3. 

We would expect by continuity arguments that the horizontal map at ^ = 1 should 
be well-approximated by following how the nearby cluster points map. Setting ra- 
tional p/q = 6 < 1, we have p ^ q, and under the LFT we have p' /q' = p/{2 *p + q). 
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Figure 9. Map from 6 = 1 line to 1/3 lines (middle interval). 



There are q points near line 1, and they map to the middle third of the 2p + q ^ 3q 
points under the LFT. 

So we do a numerical experiment: Map these cluster points near 1 to corresponding 
points near 1/3, and plot the result, to obtain Figure 9. We note the curve is almost 
linear, but not quite. Performing a best fit polynomial approximatioiiQ we find the 
map is in fact given by 



(39) 



X 



Qy-y 



We observe that these polynomial maps are known from earlier work (Choi et al. 



1990) as characteristic polynomials Pq for matrices used to compute the spectra of 



qxq matrices representing the operators hg in the rotation algebra Aq, with 6 = p/q. 
Thus, the map we observe is simply 

(40) Pi(x) = -Pi/3(y), 

where by Pp/q{x) we mean the q-th order characteristic polynomial of the qxq matrix 

(41) H = U + U* + V + V*, 

with U the cyclic permutation matrix and V a diagonal matrix with consecutive 
powers of e^'^^^^'^ on the diagonal. Details are in (Choi et al., 1990; Lamoureux and 



Mingo, 2007). 



To see the full algebraic curve, we plot the full curve in Figure 10. Observe 
the central portion of the curve corresponds to the nearly linear section we saw in 
Figure 9. 



In fact, in our first numerical experiments, we inadvertently computed the inverse of a polyno- 
mial, with expansion y — x/6 + a;^/1296 + a;^/93312 + • • • . The second author recognized this as 
the inverse of a cubic! 
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Figure 10. Algebraic curve Pi{x) = —Pi/^iy) defines the 1 to 1/3 map 



Figure 11. Map from 1/3 lines to 1/5 lines (3 middle intervals). 



8. Horizontal similarity: Quintic case 1/3 i-)- 1/5 

We repeat the numerical calculation, using the same central symmetry with LFT 
given by the matrix 



(42) 



1 

2 1 



But now we look how the three spectral lines at ^ = 1/3 maps to the three central 
spectral lines at level 9' = 1/5. A numerical experiment as in the last section produces 
the plot in Figure 11. Again, the three segments look nearly linear, but not quite. 
An inspired guess suggest these segments come from a portion of the algebraic curve 

(43) Pi/3(X) = -Pi/5(l/), 

and a comparison of the two plots shows this is indeed the case. 
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Figure 12. Algebraic curve Pi/^ix) = -Pij^iy) defines the 1/3 to 1/5 
map. Note the three nearly linear segments across the centre diagonal, 
matching Figure 11. 



It is worth noting that we can restrict the algebraic curve to the segments shown 
in Figure 12 by plotting only those points (x,2/) such that Pi/3(x) = —Pi/^{y) and 

l^i/3(a:)|<4. 

9. HORIZONTAL SIMILARITY: ALGEBRAIC CURVES 

We have seen in Sections 7 and 8 that the horizontal maps appear as nearly linear 
sections of algebraic curves, in the form 

\Pe{x) = ±Pe'{x')\<4, 

for the two characteristic polynomials Pg and Pg>, where 6,6' are rational. The full 
algebraic curves are interesting in themselves, so in this section we present a few 
plots of the curves and point out some obvious patterns. 



Referring to Figure 13, we see the algebraic curves Pg{x) = Po'iy) for rationals 
6 = 1/q \-^ 6' = l/(g + 2), with q odd, have a particularly simple graphs. There is a 
single connected component, that spirals around the origin, with more spirals as the 
denominator q increases. We see q nearly linear segments, which are the parts of the 
graph that define the relevant horizontal map on the fractal. We also see the odd 
symmetries of the curves, that point {x,y) is in the curve if and only if (— x, —y) is 
in it too. 
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Figure 13. A sequence of related algebraic curves, for rationals: 1 i— ?■ 
1/3, 1/3 ^ 1/5, 1/5 H-> 1/7, 1/7 ^ 1/9, 1/9 ^ 1/11. Notice the odd 
number of nearly linear segments down the diagonal of each map. 



Referring to Figure 14, we see the maps for rationals 1/q i-> l/(g + 2), with q even, 
have somewhat more complex graphs. There are q/2 connected components, that 
form something like figure eights around the origin, with more figure eights as the 
denominator q increases. We see q nearly linear segments, which are the parts of the 
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Figure 14. A sequence of related algebraic curves: 1/2 i— )■ 1/4, 1/4 ^^ 
1/6, 1/6 ^-)■ 1/8, 1/8 ^-)■ 1/10. Notice the symmetry, and the even 
number of nearly linear segments down the diagonal of each map. 



graph that define the relevant horizontal map on the fractal. We also see the even 
symmetries of the curves, that point {x, y) is in the curve if and only if the points 
(±a;, ±y) are in it too. 

Figure IlS] examines the curves for rational 6 of the form 6 = p/3. In both cases 
{p = 1/3,2/3), there is odd symmetry, there are three nearly linear segments, and 
some spiraling of the curve about the origin. But for the case 6 = 2/3, there are also 
two additional disconnected components to the curve, that seem to have nothing to 
do with the linear segments that are involved in the fractal map. 

In Appendix 3, we include more plots of these algebraic curves. We note with 
even denominator q, we have even symmetry in the plots, there are q nearly linear 
segments, and we see nested "figure eights" defining the key parts of the curve. For 
J9 > 1, there are disconnected components of the curve that seem to have nothing to 
do with the main part of the curve that defines the fractal map. 
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Figure 15. Algebraic curves with denominator 3: 1/3 i— ?■ 1/5, 2/3 i— ?> 
2/7. Notice the the three nearly linear segments down the diagonal of 
each map, and the disconnected components on the right graph. 

We also observe the curves appear to get much more complex with increases in the 
numerator p, for the parameter 6 = p/q. There may be some connection between 
the various straight line segments in curves of different 6 = p/q with the same 
denominator q, but these patterns are not entirely clear. 

We leave discussion of these curves and their connection with symmetries to future 
work. 



10. The SIMILARITY MAPS: GENERAL GASE 

The work in Sections 7 through 9 suggests an obvious candidate for the horizontal 
maps within spectral intervals: namely, a correspondence of points determined by 
two characteristic polynomials. Combining this observation with the earlier sections 
discussing the vertical maps and interval mapping, we can now specify the form of 
the general similarity maps of the butterfly, as follows: 

Fix a matrix element M in GL2(Z), represented as 



(44) 



M 



a b 
c d 



Fix an integer r > 0. The similarity map S 
described as 



S 



M,r,+ 



on the rational butterfly is 



(45) 



ix,e) ^Six,9) = ix',9'), 
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where the vertical component of the map is given by 

(46) e^e' "^ + ^ 



cd + d' 

while the horizontal component is given by the interval and polynomial correspon- 
dences. Specifically, recall the spectrum at level 6 = p/q is a set of q intervals Ik, 
k = l...q, while at level 6' = p'/q' there are q' intervals J^,, k' = l...q'. The map from 
X at level 6 to x' at level 6' is determined by the interval and polynomial conditions 

(47) xeh^x e /;.p,+, = h' with {-iy+''Pe{x) = i-lf+'^'Peix'). 

The sign choice for the polynomials is made so that ±Pg{x) and ±Pg'{x') are both 
monotonically increasing on Ik, I'k,, respectively. These two conditions determine the 
image point x' uniquely, except possibly in one special case. Except for this case 
(described next), we have determined the similarity maps precisely. 

The special case to consider is when x = 0,6 = p/q with q even. In this case, the 
two intervals Iq/2,Iq/2+i overlap at x = and may map to two disjoint intervals J^, 
and I'k,_^i- Here, we must decide where the point x = maps to, either an endpoint 
of I'f., or to an endpoint of I'ki^i- It is convenient to chose to map to both points, 
making the function S double- valued at this point {x = 0,6 = p/q). We will see in 
the next section that the double-valued function is continuous. 

There is another set of similarity maps, S = SM,r- , which differs from the above 
example by the way the interval maps are selected. For this maps, we chose 

(48) X e Ik^X e I'r.(q,^.p>)+k = h', 

which we observed in Section 6 as similarities that can occur. 

Of course, not all such maps SM,r,± necessarily appear as similarities of the rational 
butterfly. For instance, from Theorem 1, we know the matrix M must be chosen from 
the semigroup G C GL2(Z)/±/ described in Theorem 1. The integer r must be small 
enough that 

(49) r ■ p + q < q 
in the "-I-" case, and 

(50) r-{q'-p') + q< q' 
in the "-" case, for all 6 = p/q. Equivalently, we have 

,,^. . . c6 + d-l 

51 r < mm 

^ ^ ~ a6 + b 

in the "-I-" case, and 

(52) r < mm — — - 

e [c-a)6^d-h 
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in the "-" case. But these are the only apparent restrictions on M,r, ±. 



11. The similarity maps: Gap labelling 

The spectral gaps forming the butterfly are conveniently labeled using a Diophan- 
tine equation, where integer parameters (s, t) are flxed, and the fc-th gap in the 
spectrum = p/q is identifled using the formula 



(53) 



k = t*p~s*q. 



This indexing is known as gap labelling, as described in (Bellissard, 1990 Goldman, 



2009 Kaminker and Putnam, 2003 Ypma, 2007), where the parameters (s,t) are 



related to Chern numbers. These parameters give a convenient way of labelling the 
"wings" in the butterfly. Note there is a limited range on the indices: for t > we 
require 



(54) 

while for t < we require 

(55) 



< s < t-1, 



t < s < -1. 



It is easy to check that the similarity map 5* = SM,r,± maps labeled gaps to gaps, 
as stated in the following: 



Theorem 2. The similarity map S = SM,r,- 
{s',t') according to the formula 



maps the gap labeled (s,t) to the gap 



(56) 



{ad — be) 



a b 
c d 




s 
t 


+ 


' " 

r 



Also, the similarity map S = Sm,-, 
according to the formula 



maps the gap labeled (s, t) to the gap {s', t') 



(57) 



s 
t' 



{ad — be) 



a b 
c d 




s 
t 


+ 


r 
r 



Proof: We consider S = SM^r,+- Note the matrix M in the similarity determines 
the map from angle 6 = p/q to 6' = p' /q' according to the matrix formula 

(58) 

In the case where the determinant ad — be is one, we invert to flnd 
(59) p = {dp' — bq') and q = {—cp' + aq'). 



'p'' 


= 


a b 
c d 




P 
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The index r of the similarity tells us how intervals Ik are mapped to intervals Ik' , 
where k' = r ■ p' + k. Thus the k-th gap is mapped to gap k' = r ■ p' + k. With gap 
k = t * p — s * q, we write 

(60) k' = r ■ p + k 

(61) = r-p' + t*p — s*q 

(62) = r ■ p' + t{dp' — bq') — s{—cp' + aq) 

(63) = (cs + dt + r)p' — (as + 6t)g' 

(64) = t'*p'-s'*q', 

from which we read off s' = {as + bt) and t' = (cs + (it + r), as desired. 

The other cases (negative determinant, and S = SM,r,~) are similar. QED 

12. The similarity maps: Proof of continuity 

We show first continuity along the horizontal direction for the similarity maps 
described Section 10: 

Theorem 3. Given two rational parameters 9 = p/q, 6' = p' /q' , and a correspon- 
dence between points in pairs of intervals Ik,I'k' in the spectra Spec{hg), Spec{hg/) 
respectively, given by 

(65) xeh^x eh' = I'r.p.+k with (-l)^+''P,(x) = {-lY^^'Pe\x'), 

then the correspondence is a continuous bijection on each interval. In particular, this 
yields a continuous map from the set of points x G Spec{hg) into x' G Spec{hg/), with 
the possible exception at the point x = 0, when q is even. 

Proof: The polynomial maps {—l)'^^^Pg{x),{—iy^''Pg>{x') are both continuous 
and monotonically increasing on the intervals in question, and thus have continuous 
inverses. So on each interval, we have a continuous bijection. For q odd, the spectrum 
SpecQig) is a disjoint union of the intervals Ii...Iq, and so the collection extends to a 
continuous map into Spec{hg/). In the case where q is even, two of the intervals Iq/2 
and Iq/2+i overlap at the point a; = 0. So, except at this point, the correspondence 
extends to a continuous map on the union of intervals - that is, it is continuous at 
Spec{hg)\0. At the point a; = 0, the correspondence maps the one point to (possibly) 
two points, one at the right endpoint of /^ /_,_ /g and the other at the left endpoint 
of I'r.pi ^g/2+i- This double- valued splitting (although not a singled- valued function) 
is nevertheless continuous. QED 



Figure 16 indicates graphically how the bijection works on level 9 = 1/3 to 6 = 1/5. 
We note that the correspondences of the form 

(66) xeh^x'elk'= /;(,._p.)+fc with {-iy+''Pg{x) = i-iy'+'^'Pgix'), 
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Figure 16. Plot of cubic polynomial Pii^{x) and quintic Pi/5(x'). 
Match the zero crossing, move polynomials up and down, and observe 
a continuous correspondence of the zero crossings. This gives the con- 
tinuous correspondence between line spectra. 



also are continuous, as in the theorem. This covers the similarity maps of the form 

The numerical work suggests that piecing together these individual polynomial 
maps, at the different 9 levels, gives a continuous map on the whole butterfly. We 
now state and prove this as a theorem: 

Theorem 4. Suppose S = SM,r,± is a similarity of the rational butterfly, as described 
in Section 10. Then S is a continuous (single-valued) map, except possibly at the 
points (x = 0, 6* = p/q) with q even, where the map is double-valued and continuous. 

Proof: We consider the case 5* = SM,r,+- Fix a sequence of points (x„. On) converg- 
ing to point (x^,, 0^) and set 

(67) {x'nX) = S{Xn^en). 

By continuity of the linear fractional transformation, the sequence 0'^ converges to 
ff* = (06** + b)/{cO^ + d). To show continuity of S, we only need to verify convergence 
of the x'^. We will show every convergent subsequence of the image points (x^, 0'^) 
converges to the point (x'^, 0'^) = ^(x*, 6**), which, by compactness, establishes conti- 
nuity. 
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Take a convergent subsequence of the (x^,6'^) with hmit {x'g,0'^). Suppose in 
this subsequence, there are infinity many points where 9'^ = 6'^ (and consequently, 
Sn = &*)■ Applying Theorem 2, by continuity along the spectral lines at 6^,, we can 
conclude that x'^ converges to x'^. (Except possibly at x^ = 0, 6** = p/q, with q even. 
Here, the map S may be doubly valued, as discussed earher.) 

The other case to consider is when there are at most finitely many n with 6'^ = 6^. 
By restricting to the tail of the subsequence, and renumbering the subsequence, 
WLOG we may assume 6n 7^ 0^, for all n. 

Set the notation consistently, so in the domain we have 6'„ = Pn/qn, Q* = P*/q*, 
and points x„ G Ik„, x^ G /fe. , while in the range we have 9'^ = p'n/q'n: ^* = P'*/^!*^ ^^^ 
points x^ G /^/ , x'^ G /{,/ , x'^ G J^,/ . From the definition of S, we know k'^ = r ■p'n + k^ 
and k'^ = r ■ p'^ + k^. The value of k'^ is to be determined (it will equal k'^). 

We apply the trace tq^ to the spectral projection X(-oo,x„)(^e,J- By Lemma 7 in 
Appendix 2, we have 

(68) r,„(x(-oo,.„)(V)) = ^ + ^F{{-lf-^'^-PeSxn)l 

Hn Hn 

where the function F is an integrated density of states on ho- By continuity of the 



trace and spectral projection (see (Boca, 2001), Chap. 5), this quantity converges to 



(69) r,.(x(-oo,x,)(V)) = ^^^ + -F{{-lf'+'^'PeSx.)). 

q* q* 

Since the sequence of ratios Pn/q-n 7^ P*/q* converges to p*/q^, we have that g„ 
converges to infinity, so the terms above with l/g„ in them converge to zero. Thus, 
from the limit of the traces, we conclude 

(70) hm ^ = '^^^ + 1f((-1)'=*+^*P,.(x.)). 
n^oo g„ q^ q^ 

Similarly, by applying the trace to the sequence of spectral projections on the 
convergent sequence x'^ — )■ x'^, we find 

(71) ^.;(x(-oo,<)(v)) = ^ + ^ni-ir-'-''^Ped<)i 

converges to 

(72) redXi-o.yAheO) = ^ + ^F{{-lf^+'^'^Pedx'o)). 

y* y* 

Thus 

(73) hm ^ = ^ + \F{{-lfo^^'^Podx'o)). 
-^^^ Qn Q* Q* 



n— ^00 



26 MICHAEL P. LAMOUREUX, JAMES A. MINGO, AND SYDNEY R. PACHMANN 

With 

.-J,^^ Ki ^ r-p'^ + kn ^ r ■ {apn/qn + b) + kn/qn 

q'n q'n CPn/qn + d 

we take limits to obtain the identity 

(75) hm^ = r ■ ^: + ^_^lim^. 

qn Cfc** +d qn 

Replacing with the limiting values computed above, we obtain 
(76) 

q* q* c(7=K -r ct \ q^: q^: J 

Multiplying through by g^ and simplifying yields 

(77) K-^^ ni-^f'-^'-Pe'SO) = r ■ K + A;, - 1 + F((-l)'=*+^*Pe.(x,)). 

Now, if x^K is an interior point of the interval I^^ , this last equation has F taking 
a fractional value strictly between zero and one. Equating the integer part of the 
equation gives 

(78) k'^-l = r-p', + k,-l, 

from which we conclude k'^ = k'^, which tells us that the limit point x'^ and image 
point x'^ are in the same spectral interval /[^. 

Equating the fractional part of the equation gives 

(79) Fii-lf*^'^''Ped<)) = F((-l)'=*+«-P,.(x.)), 
which, from the polynomial correspondence between x* and x'^, yields 

(80) Fii-lf^^^'^Ped^)) = Fii-lf'-^'^'^PedO)- 

Since F, composed with polynomial (—1)''*^'^* P0r^{x) is strictly increasing on the 
interval /^^, the equality implies x'^ = x'^, as desired. 

Now, when a;^, is at an endpoint of interval Ik, we use gap labeling to show the 
hmit point x'^ must be the corresponding endpoint of interval /^^ . Suppose x* is the 
left endpoint of interval 7^: in this case, equation 77 reduces to 

(81) fco - 1 + Fii-lf''^'''Pe'Sx'o)) =r-p', + K -1 + 0. 

so we know that k'^ equals either k[ — 1 or k'^,. We just need to show it is equal to k'^. 
With the point x,,, 7^ at a left endpoint, there must be a gap to the immediate 
left of x^, which has some label (s,t). By continuity of gap labeling (Prop 11.11 
in (Boca, 2001)), the gap with this label forms an open set, so for n sufficiently 
large, the point x„ is to the right of this gap. Hence the interval Ik„ is also to the 



= 


r 


■ P'n + kn 












> 


r 


■P'n + t* 


Pn- S*qn 








= 


r 


■P'n + t* 


{dp'n - 


Hn) - 


- s * {- 


-Cp'n 


+ m'n 
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{cs + dt + 


r) * p; 


— {as 


+ bt)* 


q'n 
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t' 


*Pn-S 


*?n 
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right of the gap. Since the interval number must be bigger than the gap number, we 
have 

(82) kn> t*Pn- S*qn. 

The image point x'g, is in interval /^ and its index satisfies 

(83) K 
(84) 
(85) 
(86) 

(87) 

where the indices {s',t') label the image gap, as in Theorem 2. Thus the interval /^, 
is to the right of the gap (s', t') and hence so is the point x'„. Since this gap is open, 
the limit point x'^ is to the right of the gap, and thus the interval Ij. containing point 
x'g is to the right. Consequently, we have 

(88) K > t% - s'q: = A;: - 1. 

This eliminates the possibility that /c^ = /c^ — 1, so we are left with k'^ = k'^, from 
which we quickly conclude that the point x'^ is the left endpoint of interval /^^ and 
thus is equal to x'^. 

Handling a right endpoint is similar. A similar analysis of the case a;,^ = at the 
common endpoint of two overlapping intervals shows the image sequence x'^ could 
converge to endpoints of two different, disjoint intervals. QED 

It would seem the result on the rational butterfly should extend by continuity to 
the full butterfly, including irrational values for vertical parameter 6. We state this 
precisely as a conjecture. 

Conjecture 5. Given a similarity S = SM,r,± of the rational butterfly, as described 
in Section 10, there is a unique continuous extension of S to the full butterfly, possibly 
double-valued at certain points along the vertical line a; = 0. 

We do not have a proof of this result. It seems a density argument and use of the 
continuous trace on the fleld of rotation algebras A0 should lead to the result. 

It is curious to consider how the polynomial correspondences on the rational spec- 
tral lines might extend to irrational values of 6, in which case there are no flnite 
polynomials to describe the mapping. 
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13. Three generators for the butterfly similarities 

There are three similarity maps of the butterfly which apparently generate all the 
similarities discussed above. There is the horizontal flip H deflned as the map 

(89) H{x,e) = {-x,e). 

There is the vertical flip V deflned as the map 

(90) v{x,e) = {x,i-e), 

which can be represented in the form discussed in Section 10 as 
(91) V = Sb,o,+, 

-1 1 



with matrix B — n i 

tion 4, taking the butterfly to the bottom half given by 



for matrix A 



. Finally, there is the similarity map discussed in Sec- 

(92) S{x,e) = {x',^), 

where x G /fc in the k-th interval at level 6 is mapped to x' G /^ in the k-th interval 
at level 6'. This map is represented as 

(93) S = Sa,o,+, 

1 
1 1 

Some elementary calculations show how these three similarities combine alge- 
braically. First, there are the two obvious identities 

(94) H^ = 1, V^ = /, 
and the commutation relation 

(95) HV = VH. 

The horizontal flip almost commutes with similarity S\ in fact it just gives a shift 
in the interval indexing as follows: 

(96) HSa,o,+H = Sa,i,+. 
This shifting extends to powers of S, so we flnd 

(97) {HSAfl,+H)"' = SA",n,+ - 

By composing powers of S and HSH, we obtain all similarities of the form 

(98) 5A",r,+ for < r < n. 

We can also verify that such operators combine in the form 

(99) WA"' r' +)('-' A",r,+) — ^A"+"' r+r' +• 
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The vertical flip does not commute with similarity S, but instead introduces the 
similarity Sm^t- in the mix. Again, a straightforward calculation shows 

(100) VSAr^,T,+y = SBA-^B,r~- 

As noted in Theorem 1, the matrices A, B generate a large semigroup of linear 
fractional transformations. From the calculations noted above, the three similarities 
of the butterfly generate enough similarities to cover this semigroup. 

Theorem 6. The three similarities H, V, S of the rational butterfly generate a semi- 
group of continuous (possibly double-valued) similarities of the butterfly, including 
maps of the form 

(101) SM,r,+ and SM,r- 

where the matrices M range over the semigroup of elements of GL2{Z)/ ± I repre- 
senting linear fractional transformations mapping the interval [0,1] into itself. 

The double-valued character is, of course, from the interval splitting in the case of 
6 = p/q with q even. There may be other similarities of the butterfly. The point is, 
we can at least see this large semigroup from GL2(Z) appearing. 

14. Conclusions 

The Hofstadter butterfly, representing spectra of a continuous family of almost 
Mathieu operators, shows obvious fractal-like symmetry. By investigating these sym- 
metries numerically, we have catalogued the self-similarity maps of the butterfly using 
a semigroup of Mobius transformations in the vertical direction, indexed by elements 
in the matrix group GL2(Z). This semigroup is generated by two matrices. In the 
horizontal direction, the self-similarity maps are given by algebraic curves determined 
by characteristic polynomials. Properties of the algebraic curves are demonstrated 
in a series of plots. These algebraic curves show nearly linear segments, demon- 
strating an almost linear behaviour in the horizontal component of the self-similarity 
maps. We proved continuity of the similarity maps on the rational butterfly, possibly 
double- valued at points with horizontal parameter x = 0. The similarity maps are 
generated by exactly three continuous symmetries. We conjecture the semigroup of 
similarities on the rational butterfly extends to a family of continuous similarities on 
the full butterfly. 
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Appendix 1 

Proof of Theorem 1: We show the corresponding semigroup in GL[2, Z)/{±/} 
is generated by the two matrices A and B. Suppose the matrix 



(102) 



M 



a b 
c d 



is in the semigroup. In the case that 6 = 0, the determinant condition gives a 
±1; the condition that the LFT maps [0, 1] into [0, 1] reduces WLOG to a ■ 
and c > 0. This leave matrix M in the form 



±d 
d = 



(103) 



M 



1 
c 1 



A^ 



with c > 0. Hence the matrix M is generated by A alone. 

In the case that 6 7^ 0, WLOG (by mod-ing out by ±/), we can assume d > 
6 > as we know that maps to b/d G [0, 1] under the LFT. By the determinant 
condition ad — be = ±1, we have that the gcd of b, d is one, so we may apply the 
Euclidean algorithm to obtain a sequence of strictly positive quotients Qq, qi, q2, ■ ■ ■ 
and remainders 6 = Tq, ri, r2, . . . with 

d = qoro + ri 

b = qiTi + r2 

ri = g2^2 + ^3 



where Vk 
factor as 



gcd{b, d) 



Tk-i = qkr-k + r-k+i, 
1 and Tfc+i = 0. Applying these values to matrix M, we can 



' a b' 
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1 " 




ao ri 












c d 
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qo _ 




_ Co r-Q 
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1 ' 




' 1 " 


fll ^2 
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go _ 




. 1 ^1 . 


_ Ci ri _ 
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' 1 " 




' 1 " 




" flfc " 
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go _ 




. 1 ^1 . 




. 1 * 


c 




_ Cfc 1 _ 



for some appropriate values of ao,ai, 



, Co,Ci, 



Now, the matrix factors with the qj are generated by the A, B matrices, since 



(104) 





1 



A'^-^BA. 
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In the last factor, 



(105) 



Ofc 

Cfc 1 



we know that c^ is non-negative, else the LFT has a pole in [0,1], which is not 
allowed. By the determinant condition, we have a^ = ±1. In the case 0^ = 1, we 
have that the last factor appears as 



(106) 



1 

Ck 1 



^c,_ 



In the case where a^ = —1, we may combine the last two factors to observe 

"0 1 1 r -1 Ol _ To liri -llTc-l 1" 
Ig cl ~ lgJ[o lj[ cl 

1 1 r c-1 1 

1 g-1 J [ c 1 

= {A''-^BA){BA'). 

For g > 2, we are done: these last factors appear as generated by matrices A,B, as 
desired. 

In the case g = 1, it easy to check that we must have more than one g-type matrix 
in the factorization, for if not, we get 

(107) 

where go = 1, giving an LFT mapping 1 to (c + l)/c, which is outside the interval 

[0, 1]. With a second g-type matrix in the factorization, we get 

(108) 

M = ■ ■ ■ {A''-^BA){A-^BA){BA^'^) = ■■■ {A''-^B){AA-^){BA){BA"^) = ■■■ A''BA'\ 

which puts M in the semigroup generated by A, B. QED 

Here is a nice example, to show that the a^ = — 1 case really does appear. 

(109) 



a b 




' 1 " 




'-10" 


c d 




_ 1 go _ 




_ Co 1 _ 



2 1 

3 2 



1 

1 2 



-1 
2 1 



The matrix on the left gives an LFT which maps to 1/2, and 1 to 3/5. This simi- 
larity map actually appears in the Hofstadter butterfly, as we can see by examining 
Figure 1. We get ai = —1 in the factorization, which perhaps was not expected. 
Thus that little factor matrix on the right is not in the semigroup, since it maps 1 
to -1/3, and thus does not map interval [0, 1] to itself. 
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However, the original matrix on the left is actually in the semigroup, as shown in 
the proof. 

Appendix 2 

We prove a basic result concerning the trace of certain spectral projections in 
the rotation algebra Aq. The machinery for this result is standard, and we borrow 
heavily from the results in ( Boca, 2001[ ). 

Lemma 7. Let tq and r^ be the (framed) tracial states on rotation algebras Aq^Aq, 
respectively, with 9 = p/q. Let [a,b] be the k-th interval Ik in the spectrum of the 
operator kg = u + u* + v + v* in algebra Aq. Then for x G [a, b], 

(110) r0{x[a,x){.he)) = -ToiX[-4,i-i)'=+iPe{x))iho)), 

where X[a.x) is the spectral projection onto interval [a,x). Consequently, 

k — 1 1 

(111) Te{x(-oo,x){he)) = h -ro(x[_4,{-i)fc+<7P,{x))(/io)). 

Note: the choice of the sign in front of polynomial P0{x) ensures that the map 

(112) X ^ i-lf^'^Peix) 

is monotonic increasing on the interval Ik- In the earlier sections of this paper, it is 
convenient to write the trace-spectral projection function as the cumulative density 
F{x) for the operator ho, so we may write 

(113) Te{xi-oo,x)ihe)) = ^^ + -F((-l)'=+^P,(x)), 

q q 

where 

(114) F{X) = To{X[-4,x){ho)) 

is a continuous, strictly increasing function mapping [—4,4] onto [0, 1]. 

Proof: We do a direct calculation. The trace on Ag is obtained from a representa- 
tion of the algebra in C(T^) (8> M^(C) with generators 

(115) U = Li^Uo,V = L2^ Vo, 

where the functions Li, L2 are the coordinate maps li{zi, z^) = Zi, L2{zi, Z2) = Z2, and 
the q X q matrices Uo, K, are the usual cyclic permutation and the diagonal with 
powers of e^'^*^. The trace in this frame is just 



(116) Tg = jJ2^tr, 



<?' 
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where /i2 is the usual Haar measure on T^ and trq is the normahzed trace on q x q 



matrices. (See (Boca, 2001), page 11.) 

Fix X in 7^. The matrix- valued function h = u + u* + v + v* has exactly one 
eigenvalue in the interval Ik, which moves continuously for parameter values {zi, Z2) G 
T^. This eigenvalue contributes to the trace precisely when it lies in the interval 
[a,x]. That is, we get a contribution precisely when the value of the polynomial 
{—lY^'^ Pei^x) is smaller than z\ + 'zi^ + z\ + 'Z2'. We compute the trace in terms of 
a characteristic function on the torus T^, normalized by g, so 

(117) T0{X[a,x){he)) = - j X{{-l)''+iPg{x)<zf+zTi+z^+z^i)dfJ'2- 

With y = {—l)'''^'^Pg{x), and noting that the map {zi,Z2) t-> (2:^,2:2) has measure 
preserving inverse on the torus, we can write 

reiX[a,x)ih0)) = - X{y<zl+zTi+z^+Z2'')dfJ'2 

Q Jt2 

~ ~ X{y<zi+zT+Z2+z^) "A*2 

Q Jt2 

= ~ / XiPo(y)<zi+zT+Z2+Z2) df^2 
= -T-o(X[-4,y)(M)- 

where we use the fact that Poiy) = y and so the last integral above yields the trace 
in Ao applied to the spectral projection obtained by applying X[-4,?/] onto ho- 
using the relation y = (— l)'^+'^P6i(x), we obtain 

(118) rg{X[a,x)ihe)) = -To{X[-A,{-l)'<:+'JPg{x)){ho)), 

as desired. The trace evaluated on the full interval (— oo,x) picks up an additional 
contribution of {k — l)/q from the previous k — 1 intervals in the spectrum. QED 

It is worth noting that this trace integral can be computed explicitly, the result 
involves an integral of arccos(?//2 + cos(6')). 

Appendix 3 

We conclude with some final plots of the algebraic curves which determine the 
horizontal component of the similarity maps. Properties to notice are that the curves 
have even symmetry for denominator q even, odd symmetry for q odd, and in all 
cases there is a sequence of nearly linear segments along the diagonal line y = x. It 
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is also notably that as the numerator p increases, the algebraic curves become much 
more complicated. 




Figure 17. Algebraic curves with denominator 5: 1/5 i— ?■ 1/7, 2/5 i— ?■ 
2/9, 3/5 H-)> 3/11, 4/5 i— ?■ 4/13. In each graph, there are five nearly 
linear segments down the diagonal of the graph. All but the first graph 
have some disconnected components. 



36 



MICHAEL P. LAMOUREUX, JAMES A. MINGO, AND SYDNEY R. PACHMANN 






Figure 18. Algebraic curves with denominator 7: 1/7 i— ?> 1/9, 2/7 i— )■ 
2/11, 3/7 ^ 3/13, 4/7 ^ 4/15, 5/7 ^ 5/17, 6/7 h^ 4/19. In each 
graph, there are seven nearly linear segments down the diagonal of the 
graph. All but the first graph have some disconnected components. 
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Figure 19. Algebraic curve with denominator 2, the 1/2 f-> 1/4 map. 





Figure 20. Algebraic curves with denominator 4: 1/4 i— )■ 1/6, 3/4 i— )■ 
3/10. In each graph, there are four nearly linear segments down the 
diagonal of the graph. The right graph has some components which 
are disconnected from the linear segments. 
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Figure 21. Algebraic curves with denominator 6: 1/6 i— ?■ 1/8, 5/6 i— ?■ 
5/16. In each graph, there are six nearly linear segments down the 
diagonal of the graph. The right graph has some components which 
are disconnected from the linear segments. 
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Figure 22. Algebraic curves with denominator 8: 1/8 i— ?■ 1/10, 
3/8 ^ 3/14, 5/8 ^ 5/18, 7/8 ^ 7/22. In each graph, there are 
eight nearly linear segments down the diagonal of the graph. All but 
the first graph has some components which are disconnected from the 
linear segments. 
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